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Strong and weak coupling of two coupled qubits 
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I investigate the dynamics and power spectrum of two coupled qubits (two-level systems) under incoherent 
continuous pump and dissipation. New regimes of strong coupling are identified, that are due to additional paths 
of coherence flow in the system. Dressed states are reconstructed even in the regime of strong decoherence. The 
results are analytical and offer an exact description of strong-coupling in presence of pumping and decay in a 
nontrivial (nonlinear) system. 



I. INTRODUCTION 



Quantum information processing iQl] presupposes a coher- 
ent coupling of the fundamental bricks of quantum informa- 
tion, the qubits. In real systems, however, decoherence, dissi- 
pation and incoherent coupling to the environment is unavoid- 
able In the field of cavity quantum electrodynamics lUt], 
the notion of coherent coupling is known as strong coupling, 
in this case, between light and matter (e.g., the excited state of 
an atomic transition |4fl or an electron-hole pair in a semicon- 
ductor |5ii|). This leads to a quantum superposition of the bare 
states, resulting in so-called dressed states |@]. This regime is 
reached in systems of very high quality and under tight exper- 
imental control, so that intrinsic sources of decoherence are 
minimized as much as possible and coherent dynamics takes 
over The simplest description of strong coupling neglects dis- 
sipation altogether and thus reduces to that of mere coupling 
with strength g, introducing the notion of Rabi splitting |0]- 
Next step in the description includes the decay yi of the bare 
states, i = 1, 2 (Q]- This gives rise to a widely known criterion 
for strong-coupling: 4g > |7i — 72|- This is the case of vac- 
uum Rabi splitting where at most one excitation is involved. 
At this level, which is the most natural and fundamental since 
it describes one particle, there is no difference from the under- 
lying theoretical model. Differences appear at the next step 
of description when the excitation scheme is taken into ac- 
count. A typical description of excitation is to consider an 
initial condition. The coupling is then studied as the sponta- 
neous emission from this initial excited state. Unless the ini- 
tial condition is restricted to one excitation (as previously), the 
underlying theoretical model becomes determinant. Another 
important description of the excitation process is that of a con- 
tinuous pumping, for instance a coherent excitation that drives 
the system or an incoherent pump that feeds excita- 

tions at a given rate but without any coherent input ifTll [l2ll . 
The latter is more directly related to the intrinsic dynamics of 
the system, and is the one that will be considered in this text. 
With non-negligible pumping, the underlying theoretical de- 
scription cannot be ignored, and taking it into account leads 
to strong deviations from the paradigm of strong-coupling as 
established by the spontaneous emission of one excitation (in 
any model) lfl3ll . 
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FIG. 1 : (Color online) Schema of the systems of interest in this text. 
The main object of study, two coupled qubits, is sketched in (c). It 
will be compared throughout with other coupled systems, (a), (b), (d) 
and (e). 



An immediate extension of light-matter coupling in the 
linear regime is the linear model, namely, that of two har- 
monic oscillators with no restriction in the number of parti- 
cles Fig. [TJa)- Its comprehensive description under incoher- 
ent pumping was given in Ref. lITill . This describes for in- 
stance exciton-polaritons in planar semiconductor microcavi- 
ties is]] (where both excitons and photons are bosons). It was 
shown in this work how pumping calls for extended defini- 
tions of strong-coupling, essentially requiring absorption of 
the pumping rates in the decay rates. A more important the- 
oretical model, known as the Jaynes-Cummings model ifTsll . 
describes the coupling of a two-level system (such as an atom 
or a zero-dimensional exciton in a small quantum dot) with a 
boson mode (typically, cavity photons). Fig. [lib). It is more 
important because more closely related to a genuinely quan- 
tum regime, the linear model being essentially a classical de- 
scription cast in quantum-mechanical terms. Its description 
under incoherent pumping was given in Ref. (ill, but en- 
countered various difficulties to offer a complete picture. In 
particular, dressed modes exhibit complex patterns and a defi- 
nition of strong-coupling in this system is much more difficult 
to achieve, since splitting of the dressed states depends on the 
excitations. In particular, strong-coupling can be enforced by 
pumping, leading to situations of mixed weak and strong cou- 
pling, where some of the states are bare while some others are 
dressed. 

In this text, I will address the case of two two-level sys- 
tems (Fig.fTJcde)), which compromises between simplicity of 
the linear model and richness of the Jaynes-Cummings model. 
In particular, thanks to the reduced size of the Hilbert space, 
I will be able to solve the problem fully analytically, as in 
the linear model, a convenience not afforded by the Jaynes- 
Cummings model (when including incoherent pumping). This 
will allow me to provide a complete picture of weak and 
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strong coupling in a nontrivial system, and therefore shed Ught 
on more complicated systems. 

My description will address more particularly independent 
qubits (Fig. [Tfc)), for instance superconducting (Josephson) 
qubits I17l] . in the sense that their commutation rules will not 
be those of two fermions, that anticommute (Fig.[TJd)). I will 
also address the latter case for comparison and completeness, 
and obtain the elegant result that the expressions describing 
two coupled fermions are essentially identical to those de- 
scribing two coupled bosons, although these two systems are 
very different in character and behavior (for instance bosons 
accumulate arbitrary number of particles whereas fermions 
saturate at at most one, a distinction recovered in the formal- 
ism by merely substituting effective parameters). Also be- 
cause two two-level systems can be mapped to one four-level 
system (Fig. [lie)). I will address this case in detail, finding 
another fundamental case of interest. 



The main results, however, and the deepest connections to 
be made with other models (such as the Jaynes-Cummings) 
will be obtained from the case of two qubits. Beyond its inter- 
est for the previous reasons, the study of the coupling between 
two qubits is interesting in its own right |[T8;21 1: i t comes as 
the fundamental support of entangled states |22, 23tl, to imple- 
ment quantum gates ll24ll25ll and, in this quantum information 
processing context, the natural model to investigate decoher- 
ence i26ll27ll . In this text, the two coupled qubits will allow 
us to investigate strong coupling under incoherent pumping. 



The rest of this paper is organized as follows. In section ITU 
I introduce the model and its parameters and discuss the level 
structure. In section Hill I obtain the single-time dynamics, 
that will be shown to be the same independently of the under- 
lying model. In section HVl I obtain the power spectra, that, 
in contrast with the previous section, depend strikingly on the 
underlying model. Power spectra (photoluminescence spec- 
tra in a quantum optical context) are important because this is 
where the dressed states manifest. I will analyze in detail the 
case of two qubits (Sec. |IV Aj and contrast it with that of two 
fermions (Sec. IIV Bl ) and of a four-level system (Sec. IIV CI ). 
In section[V] I describe the strong and weak coupling regimes, 
first in the absence (Sec. IV At , and then including dV B| ) the 
incoherent continuous pump. I define new regimes of strong 
coupling, proper to the two qubits system. In section [VTl I il- 
lustrate the results of previous sections with examples of some 
interesting configurations: cases where pumping effect is op- 
timal (Sec. I VI Al l and detrimental (Sec. |VI B| i for the coher- 
ent coupling. In section IVIII based on the previous results, 
I reconstruct the dressed states, uncovering an unexpected 
manifestation of decoherence — the emergence of additional 
dressed states — to be found only in the coexistence of pump- 
ing and decay along with the coherent coupling. Finally, in 
section lVinl I give a summary of the main results as my con- 
clusions. To avoid distraction in the main text, most technical 
details appear in appendices (|B]|C]|D]i along with further ma- 
terial outside the scope of this study (lAllEli. 



II. THEORETICAL MODEL 



The Hamiltonian for two coupled qubits reads: 



(1) 



where (7i .2 are the lowering operators of the qubits, with bare 
energies Oi 2- They are hnearly coupled with strength g. The 
two modes can be detuned, by a quantity A = Oi — Oh, that 
is small enough, A ^ Oi 2, so that the rotating wave approx- 
imation is justified. The Hilbert space of the coupled system 
has dimension four, with the structure 2 (g) 2. It can be de- 
composed in three subspaces (also called manifolds, rungs, 
etc.) with a fixed number of excitations: the ground state, 
{| 0,0)}, with zero excitation, the excited state of each qubit, 
{ I 1 , 0) , I , 1 ) } , with one excitation, and the state { | 1 , 1 ) } with 
two excitations. 

An important point throughout this text is the commutation 
rules in Eq. ([T]), that are those of two distinguishable systems, 
i.e.. 



GiGi = a; a- =0, 



(7,(7 



'•=1,2,, 
(y}a,^\, / = 1,2, 
0, i^j, 



[Oi , (7y] = Oi Gj - Gj (7; = , i^j. 



(2a) 
(2b) 
(2c) 
(2d) 



Note that two operators from different systems commute. 
These commutation rules for two-level systems is most com- 
monly found in the literature of the Dicke model 1281] , that 
describes a gas of two-level systems emitting in a common 
radiation field. In the case of a fermion gas, this commu- 
tation is an approximation, that is made for the simplicity 
of the algebra and that is justified for a dilute gas by the 
fact that anticommuting operators give the same final phys- 
ical results 12911 . Indeed, when the wavefunctions of any 
two fermions is weakly-overlapping, symmetrized, antisym- 
metrized and non-symmetrized results are the same Isoll . 

In our case, the two qubits are strongly interacting, which 
sets our system apart from the Dicke model and its approxima- 
tions in many respects (see appendix lAb. For reference, I will 
analyze the antisymmetrized case of two interacting fermions 
(Fig. [TJd)) completely. However, the main object of interest 
in this text is that of two commuting qubits (like in the Dicke 
model), corresponding to the case of two distinguishable (in 
the quantum sense) qubits (Fig. [Tic)). 

The level structure of Hamiltonian ([1]) is sketched in 
Fig. 121 a) (at resonance). Such a "diamond-like" configuration 
can be mapped to a four-level system (4LS), i.e., as a single 
entity, with Hilbert space structure 1 ® 1 ® 1 ® 1 (Fig. [lie)). 
This description fits, for instance, the case of a multi-level 
atom llsill or of a single quantum dot that can host up to two 
interacting excitons (electron-hole pairs) forming a biexciton 
state 1I32I1 . On the other hand, two coupled qubits, matches 
the case of two nearby quantum dots directly coupled. The 
state with double excitation is then called an interdot biexci- 
ton state ifssll . 

The Hamiltonian Hq can be diagonalized in terms of two 



intermediate dressed states, | +) and | — ), as: 

//,) = CO I -)(- I +C0+ I +)(+ I +0)11 I 1,1)(1,1 I (3) 
with eigenfrequecies: 



«i +CO2 , 
(0± = ± ' 



(Oil =(0i+(02. 




(4a) 
(4b) 



The diagonalized level structure is sketched in Fig. |2jb). 
These are the same eigenfrequecies (0± and Rabi splitting 
(given by 2^) than for the dressed states of two coupled har- 
monic oscillators up to the first manifold ifTill . This equiva- 
lence breaks in the manifold with two excitations where the 
fermionic nature of the particles reveals and only the state 
I 1,1) is permitted, as compared to three possible states in the 
second manifold of the linear model: {|2,0),| 1,1),|0,2)}. 

The dynamics of dressed states | ±) and their spectral shape 
depend on the amount of decoherence that the dissipative and 
excitation processes induce in the system. I will consider an 
incoming flow of excitations that populate the two-level sys- 
tems at rates Pi, P2 and an outgoing flow (given by the in- 
verse lifetime) at rates 71, 72, respectively. This situation cor- 
responds to an incoherent continuous pump or injection in the 
qubit that can be varied independently from the dissipation. 
The steady state reached under the pump and decay corre- 
sponds to a statistical mixture of all possible quantum states 
and is described by a density matrix p . The master equation 
of the system has the standard Liouvillian form ifsill . with the 
corresponding Lindblad terms: 



^=^p =i[p,H] 



= 1,2 



(5a) 
(5b) 

(5c) 



=1,2 



This master equation can be exactly solved given the finite 
and small dimension of the Hilbert space (which is only four). 
Within the same formalism, we can describe the spontaneous 
emission from a general initial state by solving the equations 
for vanishing pumping. 

I will note the transitions between bare states (cf. Fig.|2a)) 

as: 



Ml =|0, \=aialo2: 
h =|0,0)(l,0|=CJi-Mi, 

U2 =1 1,0)(1,1 1= (j[(7ia2, 

/2=|0,0)(0,1 |=(J2-M2■ 



(6a) 
(6b) 

(6c) 
(6d) 



(for upper and lower transition). Note that they can all be 
written in terms of the qubit operators (in normal order) since 
Oi = M, + /,■ {i= 1,2). The Lindblad terms of pump and decay 
in Eq. (|5]) are expressed in terms of their o operators. Rewrit- 
ing them in terms of the transition operators M1.2, ^1.2 leads to 
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FIG. 2: (Color online) Energy levels for the two qubits described by 
Hamiltonian (TJ. Weak coupling (a) and strong coupling (b) in the 
absence of pump are well defined in terms of the bare and dressed 
states I ±), respectively. As shown in Sec. IV Bl the SC regime gives 
rise to new regions when pump is taken into account: (c) FSC where 
dressed states remain | ±) , an (d) SSC an MC, where dressed states 
form a new set | | 0±). Only SSC, with splitting of all the new 
dressed states, is shown. MC corresponds to the case where | I±) 
have closed (both collapsing on coi ). The thickness of the levels 
represents the uncertainty in energy due to (a,b) the decay and (c,d) 
both the pump and the decay. The arrows linking the levels due to 
pump/decay are blue for the lower (/,) and red for the upper («,) tran- 
sitions. Transitions labeled A, C (involving | +) or | in green) 
occur at frequencies determined by zi , while B and D (involving | — ) 
or I 0±), in orange) are determined by Z2- The same color code is 
used in the rest of the figures to plot the decomposition of the spectra. 



cross Lindblad terms that entangle li and mi, on the one hand, 
and I2 and U2 on the other. If the four levels did not correspond 
to two qubits but to a single entity (a 4LS), such as atomic lev- 
els or a single quantum dot levels, the Lindblad terms would 
be written directly in terms of the transition operators, without 
these cross Lindblad terms. The alternative master equation is 
discussed in appendixlC](cf. Eq. ( ICll l). 

We introduce /erm/o«ic effective broadenings: 



ri = ri+Pi, r2 = 72+^2, 

r± = — 1 — > 7± = — — 



(7a) 
(7b) 



Equation ( [Tai l is to be compared with bosonic effective 
broadenings, for boson modes a and b (such as the harmonic 
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oscillators of the linear model): 

ra = ra-Pa, tb = rb-Pb- (8) 

A tilde is being used to denote the bosonic character of the 
broadening, in the sense that the pumping strength is sub- 
tracted to the decay. Pumping leads to broadening of the line 
in the fermionic case and narrowing in the bosonic case, which 
is a spectral manifestation of Fermi and Bose statistics. 

Another convenient notation for Eqs. (|7]i is the expression 
of pumping and decay rates in terms of Fs and a new parame- 
ter r which represents the type of reservoir that the qubit is in 
contact with llssll : 

7i=Fi(l-ri), 72 = r2(l-r2), (9a) 
Pi ^ Tin, Pi^r.ri. (9b) 

All possible situations (with < r, < 1), from a medium that 
only absorbs excitation, r, = 0, to one which only provides 
them, r,- = 1, are thus included in a transparent way. For in- 
stance, a thermal bath with temperature different from zero 
corresponds to r,- < 1 /2. The effect of the medium on the ef- 
fective broadening is contained in F, . 

The power spectrum for each qubit, / = 1 , 2, is defined as; 

Si{co) ^ ((j/(«)(7i(«)) ^ —9^ r r Gf\t,t + x)e'""dtdx 
2n Jo Jo 

(10) 

where 

Gp^(f,f + T) = (cj;(f)cJ,(f + T)>, (11) 

is the first order auto correlation function. s{co) describes how 
energy is distributed and is thus of fundamental interest. In a 
quantum optical context, this can be observed directly in the 
optical emission, but for generality, I will keep the terminol- 
ogy of power spectrum. For the steady state spectrum, the 
running time t is taken at infinite values (thereby removing 
one integral). In appendix |B] we make use of the quantum 
regression formula [34] in its most general form to compute 

g['' as well as other two- and one-time correlators (second 
order correlation functions are given in appendixlE]). We shall 
focus on ; = 1 in the following, without loss of generality: 

p=A,B.C.D 

(12) 

where «i ~ J si{(o)d(0, the population of qubit 1, is 
used to normalize the expression for the spectrum (so that 
jSi{(o)d(o = iy. 




from Eqs. ^ and ( fT2] i. 



The spectrum is composed of four peaks that I label p = 
A,B,C,D, each of them with a Lorentzian (weighted by the 
coefficient Lp) and a dispersive part (weighted by Kp). The 
four resonant frequencies cOp and the associated broadenings 
(full-widths at half maximum) Yp, are intrinsic to the system, 
as they correspond to the four possible transitions in the sys- 
tem. The coefficients Kp and Lp (derived in appendix |B]) are 
the parameters that are specific to the experimental configu- 
ration (such as channel of detection) or regime (steady state 
under incoherent pumping or spontaneous emission of an ini- 
tial state). This form of S{(o) is a general feature for the power 
spectra of coupled quantum systems |3^. 

In the case of uncoupled qubits [g ~ 0), the four resonances 
reduce to the two bare energies (Oi and 0)2, broadened by the 
effective decay rates Fi and F2, respectively. The peaks are, 
in this case, pure Lorentzians. In the opposite case of very 
strong coupling {g ^ Y,P), the spectrum is also well approxi- 
mated by Lorentzians, but with the resonances cOp now at the 
dressed state frequencies o±, and broadened by the average 
rates (Fi +F2)/2. Lorentzian lineshapes correspond to the 
emission of well defined isolated modes of the system Hamil- 
tonian, weakly affected by other modes. In this case, Kp w 0. 
The dispersive contribution becomes non negligible in the in- 
termediate situations when dissipation and decoherence (or 
dephasing, cf. appendix|A| are of the order of the direct cou- 
pling. The Hamiltonian eigenmodes are then no longer neatly 
leading the dynamics and, as a consequence, their broad emis- 
sion lines overlap in energy, producing interferences. This is 
the regime of interest in our analysis, since new phenomenol- 
ogy appears for the dressed states and coupling regimes. With 
this goal in mind, I devote the next two sections to presenting 
the analytical expressions of all the quantities appearing in the 
spectrum, Eq. (fTSl) . 



III. SINGLE-TIME DYNAMICS 

We start by analyzing the relevant average quantities 
needed to compute the weights Lp and Kp in the spec- 
trum ( HjI ), that is, populations and coherences: 

« 1 = (a/ C7i ) , n2 = (aj (72) e M , (14a) 

«corr-(cTf(T2) GC, (14b) 

n 1 1 = (a/ C7i a,^ (J2) e M . (14c) 

Hi is the probability that qubit / is excited. The sum «i +«2> 
that can go up to two, is the total excitation in the system. 
«coir is the effective coherence between the qubits due to the 
direct coupling. «n is the joint probability that both qubits are 
excited. It is also the population of state | 1,1). If the qubits 
were uncoupled, we would have «ii = «i«2. 

It is important here to outline the difference between popu- 
lation of a mode (say, rii = (ct|<7i) for a qubit and «„ = (a^a) 
for an harmonic oscillator, cf. Fig.flja) and (c)), and popula- 
tion of a state (pio and pio for the state | 1,0) with | 1,0) ~ 
fl^ I 0,0) and \ 0,0), respectively). The population of the 
intermediate state | 1,0) (resp. | 0,1)) is given by «i — «ii 
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(resp. n2 — nu). This is also the probabiHty of having only 
one of the qubits excited. The population of the ground state 
is given by 1 — « i — «2 + 2« 1 1 . 

In the spontaneous emission case, ni, 112 and /icon- have the 
same solutions (depending only on the initial condition n'', 
«2 and «corr) than their counterpart for the two coupled linear 
oscillators, «„, «/, and «ab- However, the two models differ 
for nil! in the case of coupled qubits, it decays from its ini- 
tial value, nii(r) = e^'^^+'njj, whereas in the linear model, the 
population pn oscillates as a result of the exchange with the 
other states that are available in the second manifold (| 2,0) 
and I 0, 2)). Although populations of the modes have the same 
dynamics both in the two harmonic oscillators and the two 
coupled qubits, the underlying populations of their states do 
not. For instance, the decay from the initial condition | 1, 1) 
leads toni =nio + nii =Pio + Pii for the qubit, while for the 
harmonic oscillators, one has, «„ = {a^a) = Pio + P11 + 2p20- 
Since «i = (and also «2 = «&) the states are differently pop- 
ulated (pio ^ Pio, Etc.). 

In the steady state, all these mean values can be written in 
terms of effective pump and decay parameters, as in the two 
harmonic oscillators, but now following fermionic statistics 
(/=1,2): 



„ss 



aeff 



71 + 72 

ri+r2 



Qi, 



„ss 



(ni -"2), 



with the corresponding generalized Purcell rates 



Qi 



T2 Ti 
and the effective coupling strength 

„eff ^ g 



(15a) 

(15b) 
(15c) 

(16) 
(17) 



Finally, n\\ takes a simple intuitive form in the steady state, 

nfP2 



„ss 



ri+r2 

with tl 

bosonic modes {a and V), for which: 



"f^^i < „ss ss 
^ « 1 n2 ■ 



(18) 



This should be contrasted with the counterpart of « for two 



r=2nfnf- 



>nfnf. (19) 



This is an interesting manifestation of the sym- 
metry/antisymmetry of the wavefunction for two 
bosons/fermions, that is known to produce such an at- 
tractive/repulsive character for the correlators. Here we see 
that quantum (or correlated) averages (nin2) (with n, the 
number operator) are higher/smaller than classical (uncor- 
rected) averages {hi) {{12), depending on whether they are 
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FIG. 3: (Color online) Power spectrum, 5i(a)), from a qubit (thick 
solid black line) in the SC regime (71 = g and Jl = g/2) for the steady 
state under vanishing pump (Pi = 0.02g and Po = O.Olg). The spectra 
is composed of four peaks arising from the lower and upper transi- 
tions (blue and red thin lines, respectively). In dashed purple, the 
linear model spectrum for comparison. 



of a boson or fermion character, respectively. This provides 
a neat picture of bunching/antibunching from excitations of 
different modes that are otherwise of the same character. 

In the most general case, with pump and decay, before the 
steady state is reached, also the transient dynamics of the 
mean values n\, n2 and Wcorr for the coupled qubits maps to 
the coiTesponding averages «„, and n„/, of coupled har- 
monic oscillators, with only F — > F. We can conclude then, 
that the single-time dynamics of the qubit is ruled by a (half) 
Rabi frequency of the same form than in the boson case ifTill: 



R 



ITD 



^2-(F_+/A/2)2, 



(20) 



only with Fermion-like effective broadenings, Eqs. (|7]). Since, 
in contrast to the boson case where there is only one Rabi pa- 
rameter, another expression will arise in the two-time dynam- 
ics for coupled qubits, I will refer to Eq. ( |20l i as the single-time 
dynamics (half) Rabi frequency. 

We conclude this section by noting that the magnitudes 
studied up to now are independent on having two qubits, two 
identical fermions or a 4LS. The cross terms appearing in the 
master equation for the first case, due to the correlations in- 
duced by the incoherent processes, do not affect the steady 
state populations, as is shown in appendix |C] 



IV. POWER SPECTRA 

A. Two qubits 

The general expressions for the spectrum and two-time cor- 
relators of two coupled qubits admit analytic solutions at res- 
onance in the steady state of an incoherent continuous pump. 
From now on, we will refer always to this situation and, there- 
fore, I will drop the steady state label in the notation. 
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The four coefficients Lp + iKp appearing in Eq. ( fT2] i. now 
defined in the steady state, read: (cf. appendix iBb 



«l 



La + iKa {2(2^1 + i^2){R - 'T- 

+ 2^ [ - ^ (2zi + /r2) +2{R + z, + iT+)\ ^ I . 
L i + J «1 J 

(21a) 



+ '■''^B =T7^|2(222 + iT2){R + 'r_) - fli - fl2 — 

(2z2 + /r2) + 2(/?-z2-/r+)l^| 



\6Rz2 

-2g 



(21b) 



Lr + //Tr = 



^f2(2zi-/r2)(7;-/r_)-fli-fl2- 

(2zi-/r2)-2(/?-zi + /r+)l^V 



16fei 
+ 2^ 



r+ 



(21c) 



Ld + iKo {2(2z2 - iT2){R + iT- 



16Rz2 
2g 



-fll +fl2 — 



(2z2-/r2)-2(/?+z2-/r+)l^| 



(21d) 



They are defined in terms of the parameters 



ai = ^ [4r; +2Pi (P2 -2r+) -P2ri] , 02 = ^Pi {Pi - 7i ) , 

(22) 

and the corresponding frequencies and decay rates, that also 
appear explicitly in Eq. (fT2l i: 



7a 



i(OA = 2r+ + izi , 



^ + /a)c = 2r4 



Yb 



Yd 



- i(OB = 2r+ + iz2 , (23a) 



y +/C0D = 2r+-/Z2. (23b) 



They all depend on two complex parameters, z\ and Z2' 



Zl7 ■ 



-(/r+±7?)2. 



(24) 



The degree of symmetry, D^, is a real dimensionless quantity, 
between and \/2, given by 



ViYiP2 + Y2Pi)/2 



(25) 



cooperative behavior of two coupled modes, similarly to the 
superradiance of two atoms in the Dicke model or the renor- 
malization with the mean number of photons in the Jaynes- 
Cummings Model. 

The last and most important parameter appearing in the pre- 
vious expressions is a Rabi frequency for the two-time dy- 
namics, that for coupled qubits differs from its counterpart for 
single-time dynamics (cf. Eq. (|20li): 



R = Jg^-{D^gy--rl. 



(27) 



This is the true analog of the (half) Rabi frequency of the 
linear model since this value, not its single-time counterpart, 
determines strong or weak coupling (emergence of dressed 
states). At vanishing pump, the renormalized coupling G con- 
verges to g, and both R and R^^^ converge to the standard 
expression for the (half) Rabi splitting lfl4ll : 



Ro 



,2-f_ 



(28) 



The normalized power spectrum of qubit 1 follows from 
Eq. ( fTlT l with the coefficients we have obtained. The positions 
and broadenings of the four peaks are given respectively by 
the real and imaginary parts of zi and Z2- Their expressions 
remain valid in the spontaneous emission case by setting the 
pumping rates to zero. 

Figure [3] is an example of the spectrum Si{(o) (in solid 
black) and its decomposition in four peaks (thin blue and red). 
The split positions of the four peaks, which indicate the sys- 
tem is in the SC regime, are marked with two vertical blue 
lines. The two peaks that correspond to the lower manifold 
transitions, A and D in Fig. [2b), appear with a thin blue line. 
Upper transitions, B and C in Fig. |2|b), appear with a thin 
red line. All resonances are at the same positions but the 
stronger dispersive part of upper transitions leads to a shift 
of their maximum. The upper transitions are much weaker 
in intensity (magnified x30 to be visible) due to the small 
pump. The double excitation of the system is very unlikely 
(«ii = 0.0004). The lineshape is therefore close to that of 
two coupled harmonic oscillators, plotted with a dashed pur- 
ple line for comparison. The system is in the linear regime 
where all models of Fig.[T]for two coupled modes converge. In 
the following sections, we will see how the hneshapes change 
when entering the nonlinear regime. 



This quantity is proper to the coupled qubits case and its 
physical meaning will be clarified later. Its value is linked 
to the symmetry between the different parameters. For in- 
stance, Z)^ = 1 when all parameters are equal to each other, 
— Y2 — Pi = P2- On the other hand, = if one of the pa- 
rameters (any of them) is much larger than the others. It leads 
to a renormalized coupling strength 



G = D' 



(26) 



that reaches a maximum when the parameters are such that 
= V2- Such an enhancement, by V2, is related to the 



B. Two fermions 

As noted before, although the expressions for the single- 
time dynamics (populations, coherence. Etc.) for two coupled 
fermions (anticommuting operators for modes 1 and 2) are 
the same than for the two coupled qubits, their power spectra 
are different. As compared to the coupled qubits, the coupled 
fermions spectrum assumes a simple and fundamental form, 
closely related to that of the two harmonic oscillators: the for- 
mal expression is the same, differing only in the parameters 
(effective broadenings, populations. Etc.). In particular, only 
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one Rabi parameter, the single-time Rabi frequency, R , de- 
termines both the single- and two-time dynamics. The two- 
fermions spectra are thus obtained by simply substituting the 
fermionic parameters (Eq. |7| in the expression of the linear 
model I14j1 . 

This simplicity and likeliness to the linear model stems 
from the fundamental nature of the problem: two identical (in- 
distinguishable) particles coupled linearly, obeying fully their 
quantum statistics. The two coupled qubits (or the Jaynes- 
Cummings model lfl6ll ). by mixing different types of particles 
(distinguishable modes) and therefore breaking commutation 
rules, result in the more complex description and richer dy- 
namics presented in the previous section. 



C. four-level system (4LS) 

Also in the four-level system (with no cross Lindblad 
terms), the expressions for the single-time dynamics (popula- 
tions, coherence. Etc.) are the same than for the linear model, 
and here also their power spectra are different. The parame- 
ters for the 4LS spectra are of a bosonic character, cf. Eq. ([8]): 

fi=7i-A, t2 = n~P2, (29a) 
r± = , (29b) 

R=^g2-Tl. (29c) 

The relevant parameters that characterize the coupling sim- 
plify to: 




(30a) 



zi,2=R±i^+, (30b) 

recovering the conventional strong coupling criterion based on 
one parameter only, the bosonic (half) Rabi frequency R. The 
resulting spectral structure then consists of two pairs of peaks 
sitting at ±5R(^) with: 



V. STRONG AND WEAK COUPLING REGIMES 

The standard criterion for strong coupUng (SC) is based on 
the splitting at resonance of the bare states into dressed states. 
This manifests in the appearance of T-oscillations in the two- 
time correlators and a splitting of the peaks that compose their 
spectrum. 

In a naive approach to the problem of defining strong- 
coupling in a system other than the linear model, one could 
think that the condition for SC is ^ (at resonance), 

leading to the familiar inequality, g > |r_|. However, this is 
not the case whenever pump and decay are both taken into ac- 
count. Instead, one must find the condition for a splitting be- 
tween the new eigenstates, that is, the two pairs of peaks form- 
ing the spectrum. The peaks are positioned symmetrically in 
two pairs about the origin at (o,, = ±3i{zi,2) and, therefore, 

9^(zi)^0 or 5R(z2)^0 (32) 

is the mathematical condition for SC in this system. Given 
that there are two different parameters zi and Z2 on which the 
condition relies, the SCAVC distinction must be extended to 
cover new possibilities. Thus, instead of only one relevant 
parameter, r_/g, as was the case in the linear model, SC be- 
tween two qubits is determined by three parameters: 

T_/g, T+/g and D^ (33) 

This gives rise to the situations listed in Table U that are 
discussed in the following sections. 



R 




^{Z2) 


Acronym 


Type of coupling 


m 


/o 


7^0 


FSC 


First order Strong Coupling 


i\R\ 




7^0 


SSC 


Second order Strong Coupling 


i\R\ 







MC 


Mixed Coupling 


i\R\ 








WC 


Weak Coupling 



TABLE I: Type and nomenclature of coupling for two coupled 
qubits. Beyond the usual weak coupling (WC) and strong coupling 
(here denoted FSC) encountered in the linear model, the system ex- 
hibits two new regions: Mixed Coupling (coexistence of weak and 
strong coupling) and Second order Strong Coupling (with two differ- 
ent splittings of two pairs of dressed states). 



Ya 


h iCOA = 


m- 


1-/^2) H 


-71+72 


2 




4 




Yb 


^ iCOB = 


3(7H 


-72)^ 


-P1+P2 


2 




4 




Yc 


'r iCOc = 


3(7H 


-72)^ 


P1+P2 


2 




4 




Yd 


f iCOD = 


3(Pi- 


^Pi)- 


f 7i + 72 


2 




4 





iR, (31a) 

iR, (31b) 

iR, (31c) 

-iR. (3 Id) 



Note that Yp are always positive for any combination of the pa- 
rameters, in contrast with those of two bosonic modes, where 
the system can diverge. Therefore, the values of pump and 
decay rates here are not limited, always leading to a physical 
steady state. 



A. Vanishing pump and spontaneous emission 

In the case of vanishing pump, that corresponds as well 
to spontaneous emission, the standard SC and WC hold. In 
this limit, we recover the familiar expression for the half 
Rabi frequency R,R^^^ Rq- The parameters simplify to 

zia ^ V{Po±iY+? ^Po± iY+ S. 

The positions and broadenings of the four peaks are: 



Ya 
2 



- icoA = 7+ + '^0 , 



Yb 
2 



-icoB = 3y+ + iRo, (34a) 



^ + icoc = 3y+-iRQ, ^ + i(OD^Y+-iPo- (34b) 
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From here, the associated condition for SC reduces to Rq be- 
ing real, or more explicitly g > \ Y-\, as in the linear model at 
vanishing pump [see Fig.lHb)]. In SC, the two pairs of peaks 
p = A,D and p ^ B,C sit on the same frequencies although 
they have different broadenings. Excited states have shorter 
lifetime, since each excitation can decay. From Eqs. (iMt , the 
two dressed states undergo the transition into weak coupling 
(WC) simultaneously [as in Fig.[2l;a)]. In WC, Rq i\Ro\ and 
both parameters zi 2 become imaginary, giving 



C0p = 0, 
Ya 



p=A,B,C,D, 



2 =7+-\Ro\, ^ = 3r+-\Ro\, 



2 



3r+ 



\Ri, 



Yd 
2 



(35a) 
(35b) 

(35c) 



with Yp > 0, since 7+ > \Ro\ in this regime. The four peaks 
collapse into four Lorentzians at the origin, all differing in 
their broadenings. 

As a result of the two pairs of peaks sitting always on the 
same two (or one) frequencies, the final spectra can only be 
either a single peak or a doublet, both shapes being possi- 
ble in SC or WC regimes (as in the linear model and for 
the same reasons ifTill ). An intuitive derivation and interpre- 
tation of these results is given in appendix |Dl based on the 
so-called manifold picture |[l3ll . which consists in considering 
transitions between eigenstates of a non-hermitian Hamilto- 
nian, with energies broadened by the imaginary part. 

In the steady state case but in the limit of vanishing pump 
(the linear regime), only the vacuum and first manifold are 
populated. The spectra in this limit converge with the lin- 
ear model and also it can be analyzed in terms of manifolds 
by straightforward extension. The spectrum in Fig. (|3]l is an 
example of SC for vanishing pump as we can see from the 
fact that the lower transition peaks, in blue, dominate over the 
broader and weak upper peaks, in red. In this case, the split- 
ting of the dressed modes gives rise to a splitting in the final 
spectrum (in black). In what follows, we take this SC config- 
uration (71 = g and 72 = g/2) as a starting point to explore the 
effect of a non-negligible incoherent continuous pump. 



B. Non-negligible pump 

When pump is taken into account, all the types of coupling 
listed in TableUare accessible. These are plotted in Fig.|4]as 
a function of the pumping rates. The starting point is the Rabi 
frequency R, Eq. ( |27] i. that is either real or pure imaginary. 

Fkst, let us consider the case where: 



R^\R\ 



(36) 



from which follows that zi =29, and therefore 5R(zi ) = 5R(z2)- 
This is the most standard situation that we already found in 
the absence of pumping. It is sketched in Fig. |2jc). To dis- 
tinguish it from the other types of coupling to be discussed 
shortly, I will from now on call it First order Strong Coupling 
(FSC). Note that condition (|36]l can only be satisfied if < 1, 
therefore, when the renormalization of the coupling through 



5 

CM 



1 




--g/2 


1 


wcv 




MC 


\ FSC 

1 








wc; 


\ 


r MC 









2 3 

Pi/g 



FIG. 4: (Color online) Phase space of the steady state StrongAVeak 
Coupling regimes as a function of pump for Y\= g and 72 = g/2. In 
Strong Coupling (SC, blue), one can distinguish two regions. First 
order (FSC, light blue) and Second order (SSC, dark blue) Strong 
Coupling. Weak Coupling (WC) is in purple and Mixed Coupling 
(MC) in green. The dashed blue lines enclose the two regions where 
two peaks can be resolved in the power spectrum of the first qubit, 
S\ (co). One falls in SC and the other in WC. 



the interplay of pump and decay is detrimental (since Eq i 
implies that G < g). It is not possible to reach the optimum 
effective coupling and maximum splitting of the spectral lines 
given by 2\/2g. Eq. (|36] | leads to the following explicit con- 
dition for g: 



> 



v/|l-WI 



(37) 



IfR^O, then, also 9?(zi) = 9^(z2) 7^ 0. FSC includes the 
standard SC regime in the absence of pump {g > |r_ | that is 
implied by Eq. (|37])). It is the most extended region in Fig.Hl 
colored in light blue. In this case, the spectrum of emission 
follows the expected pattern: two pairs of peaks. A, D and 
B, C, are placed one on top of each other, although they are 
differently broadened [see the spectra in Fig.[3l. 

In Fig.|5ja) and (b) we track the broadenings and positions 
of the four peaks (A and D in blue and B and C in red) as a 
function of pump, through the SC region of FigH] on a diago- 
nal line defined by P2 ~ Pi/2. Following them from vanishing 
pump, where the manifold picture is exact, the four peaks can 
be easily associated with the lower and upper transitions of 
Fig. |2jb), and that is why we keep the same color code and 
notation. Dressed states | ±), close to the Hamiltonian ones, 
can still be defined in the system, but with the modified fre- 
quencies (Oi ± 3i{zi ), both affected equally by decoherence. 

By construction, the resulting spectra in this regime can 
only be a doublet or a single peak, depending on the mag- 
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nitude of the broadening of the peaks (that always increases 
with pump and decay) against the spHtting of the Hnes (that 
always decreases). As in the limit of vanishing pump, ob- 
serving a doublet in the spectra does not imply splitting of the 
dressed states (and thus, SC) IstIi . but here the tendency is al- 
ways the same: the lower the pump and the decay, the better 
the resolution of the splitting. 

Second, let us consider the situation of the Rabi frequency 
being imaginary: 

R^i\R\ ^ G->g^-rl. (38) 

This results in three possibilities, listed in Table |I](WC, SSC 
and MC), that constitute the three remaining regions delimited 
in Fig. |4] In what follows we find the specificities of each of 
these three regimes. 

The Weak Coupling regime (WC, in purple) is characterized 

by 

zi=i\zi\, Z2 = i\z2\, zi^Z2 ^ G<|r+-|/?|| (39) 

and therefore = S^fe) = 0. Note that condition ( |39] l 

is not analytical in terms of the relevant parameters i33[ . In 
WC, the four peaks are placed at the origin with four different 
broadenings. The dressed states have collapsed in energy to 
©1. 

Up to here, we have remained within the SC and WC re- 
gions already known from the linear model. We now consider 
the two new regions of SC, proper to the coupled qubits, that 
I call SSC and MC, respectively: 

SSC: When both parameters zi.2 are real, then: 

zi = ki|, Z2 = \Z2\, zi<Z2 ^ G>|r+ + |/?||. (40) 

We refer to it as Second order Strong Coupling regime (SSC, 
colored in dark blue in the phase space). Here, the broad- 
enings of the four peaks are equal, 7^/2 = 2r+, but the po- 
sitions of the pairs of peaks are different, (Oa.c = ikil and 
(Osfi = ±|z2|- The reason is that the bare energies of the 
modes undergo a second order anticrossing induced by the 
interplay between coupling, pump and decay. The energies of 
the dressed states are affected differently by decoherence, up 
to the point where we may picture the physics in terms of a 
new type of eigenstates. The association of the A and D (with 
COa d) as the peaks corresponding to lower transitions and B 
and C (with (Ob.c) to upper transitions is completely arbitrary 
in this region, given that the broadenings of the peaks, which 
led us to such association in FSC, are now equal. This implies 
that, rather than two dressed states, | — ) and | +) (Fig.|2|b)) 
as in the conventional strong coupling (FSC), the system now 
exhibits four dressed states: | /_), | /+), | (9_) and | 0+). They 
are plotted in Fig.|2|d): 1 1±) (resp. | 0±)) have energies split 
at ±|zi |, giving rise to the inner peaks, in green (resp. ±|z2|, 
giving rise to the outer peaks, in orange). The physical ori- 
gin of this remarkable departure from the conventional strong 
coupling picture will be discussed in section lVTIl 

We can see how peak broadenings and positions change 
when going from FSC to SSC in Fig. |5jc) and (d). In this 
case, we track the peaks by varying P2 for a fixed Pi, mov- 
ing upwards in the phase space. The first vertical guideline 
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P2 = Pi/2 (a) 




Pi = 0.2fl ^^^-^ (C) 
^ 


^^^^^^^^^^^^^^^ 






^ ■ (b) 

FSC 

\ ^ — — — WC 




\ MC 









01234 50 12345 
Pi 1 9 Pits 



FIG. 5: (Color online) Broadenings (a), (c), and positions (b), (d) of 
the lines that compose the spectra as a function of pump for the decay 
parameters yi = g and fi = In the plots of the first column, 
the pump P\ varies with P2 = P\/2, moving upwards in the phase 
space of FiglD The vertical guideline shows the crossing from FSC 
to WC. In the plots of the second column, the pump P2 varies with 
P\ = 0.2g, moving in diagonal in the phase space of Fig |4] The 
vertical guidelines show the crossing from FSC to SSC and finally 
to MC. The dashed blue line represents the splitting as it is resolved 
in the final spectrum Si (ft)). The color code (red-blue and orange- 
green) corresponds to that of the transitions in Fig. [2] 

marks the border between the two kinds of SC, with the open- 
ing of a "bubble" for the positions 0^ and cob (that were equal 
in the FSC region), and the convergence of all the broaden- 
ings. In principle, one can expect that quadruplets and triplets 
may form out of the four peaks. However, the broadenings 
and contributions of the dispersive parts (given by Kp) are too 
large to let any fine splitting emerge clearly. The spectra in 
this region reduce to singlets and doublets. However, we show 
in Sec.[yi]through some examples that they may be distorted, 
doubtlessly reflecting the multiplet structure. 

MC: When zi is imaginary and Z2 real, or equivalently, 

zi^i\zi\, Z2^\Z2\ ^ |r+-|/?||<G<|r+ + |/?||, 

(41) 

we enter the last new region in Fig.|4] This is a Mixed Cou- 
pling regime (MC, colored in green in the phase space) where 
the two inner peaks, A and C — as well as the reconstructed 
eigenstate | I±) — have collapsed at the origin, like in WC. 
However, the two outer peaks, B and D — as well as | 0±) — 
are still split. As in SSC, the broadening of the peaks does 
not allow for a distinction between upper of lower resonances. 
The collapsed resonance is in this case at the bare energy 
Oi = because that is the total average bare resonance in the 
system. In Sec. IVI A| we will see that when the qubits are 
detuned, this resonance happens at (oi + ft)2)/2 = —A/2. 

Again, although one may expect a triplet in MC, only dis- 
torted singlets are observed in the best of cases due to the 
broadening and dispersive parts. In Fig. |5jc) and (d) we can 
see the transition from SSC into MC, at the second vertical 
line. 

Note that, in this system, the pumping mechanism is equiv- 
alent to an upward decay, due to the ultimate saturation of 
the qubit and the symmetry in the schema of levels that they 
form. The master equation is symmetrical under exchange of 
the pump and the decay (7 <-> P,) when the two-levels of both 
qubit are inverted (| 0,0) ^| 1,1) and | 1,0) o| 0,1)) ll40ll . 
Consequently, the parameters zi, Z2 and R, and also the popu- 
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lations of all the levels, are symmetric in the same way, as it 
happens with just one qubit. In other systems, like the linear 
model, the Jayne-Cummings model or simply a single har- 
monic oscillator, the effect of the pump extends upwards to an 
infinite number of manifolds while the decay cannot bring the 
system lower than the ground state. There is no natural trun- 
cation for the pump (that ultimately leads to a divergence), as 
there is for the decay. But with coupled qubit, state | 1,1) 
is the upper counterpart of | 0,0), undergoing a saturation. 
This implies, for instance, that the (transient) dynamics in the 
limit of vanishing decay is exactly the same as that of vanish- 
ing pump and that in such case we can also apply the mani- 
fold method to obtain the right positions and broadenings as a 
function of pump, in the same way that we did as a function of 
decay only. We only have to take into account the mentioned 
symmetry consistently. As long as the dynamics moves up- 
wards or downwards only, even when intermediate states are 
coupled, the manifold picture is suitable. The manifold di- 
agonalization breaks, however, in the presence of both non- 
negligible pump and decay. This is discussed in appendixlD] 

VI. PARTICULAR CASES 

In this section, we illustrate the rather abstract previous 
discussions with examples. The symmetry in the decay and 
pumping rates determines the effective coupling, emission 
properties and dressed states. Let us start by expressing D\ 
the magnitude quantifying such symmetry, in terms of the 
reservoir parameters (F, and r,), 

= Vl ^n+n-lnn . (42) 

In this form, its physical meaning is more clear There are two 
separate factors to discuss: the symmetry in the strength of 
the couplings to the reservoirs, given by \f^Wil\^\ +r2)/2] 
and plotted in Fig.|6ta), and the symmetry in the nature of the 
reservoirs given by \Jr\ + r2 — 2rir2, plotted in Fig.HJb). 
The coupling g is enhanced when 

g<G< Vlg, that is, 1<D' <V2, (43) 

which happens when ri > 1/2 and r2 < 1/2 (or the other way 
around). This corresponds to the two squared regions with 
lighter colors in Fig. |6tb). Then, the two reservoirs are of 
opposite natures: the reservoir of the first qubit provides ex- 
citations (Pi > 7i) while the other absorbs them (P2 < 72). If 
this is accompanied by similar interaction strengths, Fi ~ F2, 
enhancement occurs. We refer to these situations as optimally 
pumped and study them in Sec. I VI Al 

On the other hand, if the reservoirs are of the similar na- 
tures, both , r2 > 1/2 or < 1/2, both providing or absorbing 
particles, then the system is detrimentally pumped: 

0<G<g, that is, < < 1 . (44) 

It corresponds to the two squared regions with darker colors 
in Fig.itb). We study this in Sec. lVTB] 




FIG. 6: Factors contributing to D": (a) yT\T\/[{ri +r2)/2] as a 
function of ri,r2 and (b) + ''2 " 2?'i ''2 as a function of ri,r2. 
The values corresponding to the contour lines are marked on the 
plots. Both functions take values from (dark blue) to 1 (light). 

A. Optimally pumped cases: g <G< \flg 

Let us explore the optimally pumped cases by consider- 
ing parameters on the diagonal Fi = F2 in Fig. |6la) together 
with the antidiagonal ri + r2 = 1 in Fig.|6jb). The reservoirs 
have opposite nature but interact with equal strength with the 
qubits. This corresponds to the situation where the decay and 
pumping parameters are equal in a crossed way: 

Pi = 72, and P2 = 7i. (45) 

The system has a total input that is equal to the total output, 
Ptot =P\+P2 = 7rOT = 7i + 72, and also equal Purcell rates, 
Q\ = Qi- The excited and ground states are formally equiva- 
lent in the dynamics. 

Figure [7] shows the different coupling regimes accessible 
with this configuration, as a function of Pj and P2 with the 
same color code than in Fig.|4] This configuration is in FSC 
only when all parameters are equal, ~ \, (blue line) and 
there is total symmetry in the system. Otherwise, one of the 
new type of coupling (SSC in blue or MC in green) is realized 
as the coupling is effectively improved, G > ^. In the inset the 
type of spectral shapes that results is shown. 



Lineshape 


Asl 


^con 


singlet 


1 


2 


distorted singlet 


1 


6 


doublet 


3 


4 


distorted doublet 


3 


8 


triplet 


5 


6 


quadruplet 


7 


8 



TABLE II: The lineshapes 5i(a)) are defined by two quantities: 
is the number of times that S\{(0) changes slope, that is, the number 
of real solutions to the equation dS\ {co)/d(0 = 0; Leon is the number 
of times that Si (co) changes concavity, that is, the number of real 
solutions to the equation d-Si (co)/dar = 0. 

The vertical axis in Fig.|7] with Pi = 72 = 7 and P2 = 7i = 0, 
is illustrative of all the possible coupling regions and line- 
shapes. This is the extreme situation of optimal pumping 
where the two reservoirs interact equally strongly with the 
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P2/g = n/g 



FIG. 7: (Color online) Phase space of FSC/SSC/MCAVC as function 
of Pl/s = 72/ 8 ™d P2/g = Yi/g- The color code is that of Fig. [4] In 
inset, the possible lineshapes of 5i(a)): a doublet (red), a distorted 
doublet (green), a distorted singlet (blue) and a singlet (white), as 
explained in Table Bll 



two qubits (Fi =T2 = j) ^nd have completely opposite na- 
tures; one only providing particles (rj = 1) and the other only 
absorbing them (r2 — 0). At this point, there is maximum 
renormalization of the coupling, G = Vlg (D** = \/2) as both 
factors in Eq. (l42T i are maximum. The populations and mean 
values read: 



«2 



4+ (7/^)2 
1 

y/g 



, ni = 1-112, 



4+ (7/^)2' 



(46a) 
(46b) 
(46c) 



The two qubits share one excitation only. The Rabi frequency 
also simplifies lo R = ig (as r_ = 0), and 



21,2 



g2_(y/2)2zpgy. 



(47) 



In Fig.|8] we can see some of these magnitudes varying in the 
different regimes as a function of y/g. In the linear regime 
limit, 7 <C g, there is FSC with all the levels equally popu- 
lated («i = «2 = 1/2, nil 1/4) and «corr = -/(7/g)/4. As 
soon as there is pumping, the SSC opens a "bubble" in the 
eigenenergies with the splitting of inner and outer peaks. The 
transition into MC, with the collapse of the inner peaks, takes 
place at 7= 2(V2— l)g, and the transition into WC, closing 
the bubble, takes place at 7 = 2(V2 + \ )g. The maximum of 
Z2 = a/2§ (in orange) takes place at 7 = 2§, when the coher- 
ence jncorrl = 1/4 is maximum. This is a special point where 
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FIG. 8: (Color online) (a) Positions of the peaks (0)4 in thick green, 
0)b in thick orange) and populations {n\ in dashed blue, «2 in dashed 
purple, |ncorr| in dotted brown) as function oi P\/ g = 72/ g = Y/g, for 
P2 = 71 = 0. The system goes from FSC (at 0) to SSC, to MC, to WC, 
while the lineshape of the spectra changes as coded in Table HH The 
most interesting lineshapes, that can only appear in SSC and MC, are 
the distorted doublet (I) and singlet (II). The total spectra (in black) is 
decomposed in inner (green) and outer (orange) peaks coming from 
the transitions sketched in Fig.[2td). 



the splitting of the dressed modes is the largest possible, 2^/2g 
(even though the lineshape remains a singlet). Finally, when 
the coupling becomes very weak, 7 ^ ^, the first dot satu- 
rates, quenching the exchange of excitation between the qubit 

and«2 =«11 =«corr = 0. 

Although the underlying physics of coupling is very rich 
and complex, the spectra do not acquire distinctively marked 
lineshapes (such as well resolved triplets of quadruplets): a 
doublet in SSC only gets distorted, but this is unambigu- 
ously due to the underlying quadruplet structure, as shown in 
Fig.Ol), and also a singlet gets distorted due to the underly- 
ing triplet structure, in (II). Before reaching WC, the spectrum 
has become a plain singlet. The way to distinguish mathemat- 
ically the different possible shapes and their origin in under- 
lying triplets and quadruplets is by counting zeros of first and 
second derivatives, as given in Table HH of which only the four 
first lines are realized in the coupled qubit. 

In Fig. |9]we make the comparison between the results dis- 
cussed in Fig. [8] and those obtained for a 4LS (without cross 
Lindblad terms), where the system undergoes a standard tran- 
sition SC-WC as 7 increases. From Sec. IIVCI we know 
that, although the populations do not change from those in 
Eqs. ( |46l ), the position of the four peaks differ, as well as the 
coupling regimes. In this case the parameters z converge to: 



Zl7 



-(7/2)2. (48) 
This corresponds to a splitting intermediate between those of 
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FIG. 9: (Color online) Comparison between the results of Fig. [8] 
(solid lines) and those obtained for a 4LS (dashed black lines). The 
positions of the peaks composing the spectra, are marked with ver- 
tical lines in (I) and (II). In a 4LS, only plain doublets and singlets 
arise. The populations are the same for both cases. Different types 
of splitting lead to dramatic differences in the spectral shapes. 



X^^^^^^ = g/2 




7 = .9 "X, 


7 = 2(! 




(11) 
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FIG. 10: (Color online) (a) Positions of the peaks {co^.c in thick 
green, (Ob.d in thick orange) as a function of detuning for six cases in 
Fig.[8]with the specified parameter y. The positions in the absence of 
crossed terms appear in dashed black and in the absence of coupling 
in thin black, for comparison. The first two plots are those labeled 
(I) and (II) in previous Figures. 



the inner and outer peaks, as we can see in dashed black in 
Fig.|9|a). Although the 4LS splitting is smaller than the spUt- 
ting between outer peaks (in orange), the SC doublet in the 
spectrum is better resolved than in the case of two qubit, due 
to the large intensity of the inner peaks. 

Fig. [To] displays the behaviors of all the resonances as a 
function of detuning for six cases in Fig. [8] (see insets), in- 
cluding (I) y — g/2 and (II) J = g- Inner and outer peaks 




FIG. 1 1 : (Color online) Spectra for varying positive detuning (an- 
ticrossings) in SSC (I) and in MC (II) for the corresponding cases 
in Fig. [21 The inset in (II) shows the positions of the three/four 
peaks composing the spectra as a function of detuning. Close to res- 
onance, the inner peaks converge to the average transition energy 
ft)li/2 = A/2. At large detunings, the first dot (that is pumped) emits 
at (Bi = and dominates over the second dot (that decays) with emis- 
sion at 0)2 = —A. 



superimpose whenever there are only one or two resonances. 
Otherwise, green lines represent inner resonances and orange 
lines the outer ones. Dashed black lines represent the case of a 
4LS, which simply splits in two lines close to resonance in SC 
(the transition to WC happens at 7 = 2^). In thin black, the 
bare energies (cOi = 0, 0)2 = —A) are plotted as a reference, 
and recovered in all cases very far from resonance or well 
into WC. The 4LS resonances converge to the bare energies 
at all detunings almost as soon as they enter WC. However, 
the resonances of two qubits, when closing in WC (MC in 
the case of inner resonances), collapse to the average between 
the bare energies (a)i + fi)2)/2 = —A/2 instead. At resonance 
these two behaviors are equivalent, but at small detunings they 
are clearly very different. We will discuss the reason for this 
puzzling saturation at —A/2 instead of —A in Sec. I VIII Even- 
tually, at very large 7, the bare energies will be recovered both 
by inner and outer peaks. 

The anticrossing that the lineshapes form when detuning 
between the modes is varied from zero to some detuning Amax, 
is also peculiar. In Fig. fTTT l) and (II) we can see that the dis- 
torted doublet and singlet keep their features up to A^ax = g 
and Amax = 2g (resp.). The resonances in Fig. [TT| are not 
equally present in the emission. In both cases, at large de- 
tuning, the emission at Oi = of the first qubit is dominant 
over that of the second qubit at 02 = —A because the first 
qubit is being pumped and the second only dissipates the ex- 
citation. Note that, the singlet in (II) slightly oscillates to the 
left before joining the origin at large detuning. This is a clear 
signature that the central peak is not simply a feature of WC, 
but of a more complicated MC eigenstate structure, as can be 
seen in the inset. This leftwards shift is a consequence of the 
splittings saturating at —A/2 instead of the bare energies (0 
and —A). This shift is even more evident when both inner and 
outer peaks saturate at some detuning, like in case 7 = Ag. 
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B. Detrimentally pumped cases: < G < g 

To explore the detrimentally pumped cases, we consider pa- 
rameters on the diagonal ri = r2 = r in Fig. IHb), where the 
reservoirs have the same nature. 

If the reservoirs also have the same interaction strength with 
the qubit, Fj = F2, the two qubits become indistinguishable. 
This corresponds to equal dissipation and pump. 



7i = 72 = 7 and Pj = P2 



(49) 



The symmetry in the system is not total as the income- 
outcome flows of excitation may not be balanced, Ttot 7^ 
PrOT- The steady state populations are the same as for the un- 
coupled case with «i = no ~ P/ (y+P), "ii = «i«2, «corr = 0, 
as happens with any two indistinguishable coupled modes 
(like, for instance, harmonic oscillators). This does not mean 
that the qubit are uncoupled, in fact, the system can be con- 
sidered always in strong coupling (FSC) with 



yP 



{y+P)/2 



--2^ril-r). 



(50a) 
(50b) 



If P ^ 7 (r = 1), or the other way around (r = 0), the symme- 
try between the flows (inwards and outwards) is completely 
broken and = 0, although, at the same time, the splitting in 
dressed modes is maximum R = g (although not enhanced). 

On the other hand, if P = 7 (r = 1 /2), the symmetry is total 
(D^ = 1). Here, all the levels are equally populated (ni =n2 = 
1/2, «ii = 1/4) and the Rabi frequency vanishes R = 0. The 
SC/WC regimes become conventional with zi,2 ~ — y^. 
The SC condition is simply given hy g > y. This particular 
situation is depicted in Fig.|2]with a thick dark blue line in SC 
that goes into WC when the parameters equal 1 . 

A second possibility that results in degrading the effective 
coupling, still considering reservoirs of the same nature (ri = 
r2 ~ r), is when r ~ 1/2, which means: 



7i=Pi=Fi/2 and 72=P2 = F2/2. 



(51) 



Although the parameters may be different, Fi 7^ F2, there 
is compensation between the flows (inwards and outwards) 
TroT = ^OT- This configuration arises when both qubits 
are in a thermal bath of infinite temperature. The popula- 
tions are also those of the uncoupled system, «i = n2 = 1/2, 
nil = nin2, neon- = 0. The Rabi and degree of symmetry read 



R 



g^ - Tl and D' = 



(52) 



F+ V° -- (Fi+F2)/2 

The FSC condition reads now g > F+ . 

VII. DRESSED STATES 

Now that we have described the new regimes appearing in 
this system due to the incoherent pumping, SSC and MC, and 
their spectral features, we can discuss their origin. 



We have pointed out that the manifold picture, based on the 
Hamiltonian dynamics, breaks in presence of non-negligible 
pumping. To define dressed states in such an essentially non- 
Hamiltonian system, one must recourse to another formalism, 
presented in this section, based on studying the flow of coher- 
ence in the system. This allows a reconstruction of dressed 
states that recovers the manifold picture for vanishing pump 
and extends it otherwise. 

To analyze how coherence flows in the system, let us ad- 
dress the dynamics of coherence between states, given by off- 
diagonal elements of the density matrix p . We thus consider 
the vector of the relevant transitions: 



^coh ■ 



/(«l)\ 




(/l> 1 




(«2) 




V(/2>y 





1,1 Ip |0,1) 
l,0|p |0,0) 
l,l|p|l,0) 

,(0,1 Ip |o,o), 



The dynamics of this vector is ruled by a matrix Mcoh: 

d\coh/dt = -McohVcoh 

that reads: 



(53) 



(54) 



/ 7. 



M, 



coh — 



72 



V 



-72 

-'g 





-Pi 

^+P2 





-'g 


7i - 'A 

-71 








>g 




-Pi 




^Pl- 





(55) 

The diagonal terms in Mcoh give the decay rate of each coher- 
ence and the off-diagonal terms, the rates at which coherences 
is transferred from one transition to another Its eigenvalues 
recover the positions and broadenings of the spectrum peaks 
in Eq. ( |23] |. with eigenstates that I note (Tp), p = A,B,C,D. 
Therefore, the single time dynamics of (Tp) is free: 

(Tp(f )) = e-('%+r./2)('-'o) (Tp(fo)) , (56) 

and their spectral shape is a Lorentzian: 



(7p/2)2 + (ffl-fflp)2- 

(57) 

In this way, one can reconstruct the dressed states from the 
above expressions by inspection of the (Tp). 

First, let us consider the simplest case of Hamiltonian dy- 
namics (without pumping and decay), for which we already 
obtained the dressed states through the manifold picture (by 
diagonalization of the Hamiltonian with only direct coupling 
g). The following discussion will thus reconstruct Fig. [2jb). 
The eigenenergy (Oa = (Ob = +ig corresponds to the two 
eigenstates: 

(Ta> = {di +d2) (lower transition | 0,0) o| +)), (58a) 
(Jb) = («2-«i) (upper transition | -) o| 1,1)). (58b) 

The eigenenergy (Oc = (Od = —ig corresponds to the counter- 
part eigenstates: 

(Tc) = {di - d2) (transition | 0,0) o| -)) , (59a) 
(Td) = (m2 + M 1 ) (transition | +) o | 1 , 1 ) ) . (59b) 
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Therefore, the system resonances 
of the dressed states: 



Tp) give away the existence 



±)oc±| l,0)+|0,l) 



(60) 



(left not normalized for convenience here and further-on). We 
can clearly see by inspection of Mcoh how the direct coupling 
transfers coherence between the two lower transitions on 
the one hand, and the two upper transitions m,-, on the other 
hand. In this case, coherence can follow a closed path forming 
loops that flow always at the same rate 2g (in between mi <-> M2 
and d\ o d2) and give rise to Rabi oscillations and dressed 
states. 

In the case of a 4LS, the situation is not much different, 
even when adding pump and decay: dressed states are also 
identified through the same procedure, diagonalizing the cor- 
responding matrix without incoherent crossed terms: 



/ 2 



M, 



coh 



V 



2 ' ri 











12 
2 



-Ig 



■7i- 








/A 



tg 




(61) 

In this case, the eigenenergies we found in Eq. dSTl ) corre- 
spond to the same combination of transitions, Eqs. ( |58]|59] l 
with: 



r_ ± iR 



1,0) 



0,1) 



(62) 



If f _ = (^ = g), we recover the dressed states without de- 
coherence. Up to the transition into WC at |f_| ^ g, both 
coefficients in Eq. (|62] | are complex numbers with unit norm 
and, therefore, they give rise to a balanced contribution of the 
two states | 1,0) and | 0, 1) equal to 1/2. Their mean energy, 
(± I // I ±) = is exactly the one appearing in the spectral 
resonances cOp. The straightforward identification of dressed 
states, in the same way as when there is only direct coupling, 
makes it possible to describe the dynamics in the manifold 
picture as explained in appendix iDl 

When the system is in WC {R = i\R\), the dressed states 
become essentially the bare ones, | 1,0) and | 0,1), since 
\r^\/g ^ That is, in this regime, the populations become 
less and less balanced, as we can see in the example of Fig. [12] 
where the probability to find the system in | 1 , 0) is plotted in 
thin dashed for both | ±). 

On the other hand, in the case of two qubits which are only 
pumped and decaying, with no direct coupling, the eigenstates 
of Mcoh are combinations of the transitions corresponding to 
each qubit: (P2M1 + Yidi), {P\U2 + 71^/2), («i — d\) and {112 — 
c/2) . No new dressed states arise although coherence is being 
transferred (incoherently) between two transitions. 

From these three limiting situations, we learn that direct 
coupling is essential for the appearance of dressed states but 
also that pump and decay can transfer coherence in the direc- 
tions Ml -H> d\ and M2 -H> d2, in an incoherent continuous way. 
In spite of this, they do not induce oscillations in the dynam- 
ics. Even in the most symmetrical case where 71 — Pi and 
72 — P2, pump and decay rather move around coherence in an 
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FIG. 12: (Color online) (a) Contributions from state | 1,0) to the 
dressed states appearing in an optimally pumped case: P\ = Ji = y, 
P2 = y\= 0. The X-axis is y/g, as in previous Figs.[8]and [9]describing 
the same case. In thin dashed black lines, the contribution to the 4LS 
dressed states | ±) . In thick lines, the qubit dressed states: orange for 
I 0±) and green for | /±). The splitting of each pair of lines ± from 
equal to different contributions, marks their transition into WC. All 
the "— " states tend towards | 1 , 0) in the very weakly coupled regime. 



|1."> 10.1) 

h 



|1,0) 10,1) 



9 . 



l2\ 





"1 




|1,0) 




10.1) 



< — ► 



|1,0)1 10,1) 



FIG. 13: (Color online) Circulation of coherence among the transi- 
tions between bare states for two qubits in the case where y\=P2 = 0. 
A loop is formed that involves all the transitions, thanks to the pump- 
ing and decay terms that participate to the flow of coherence in a con- 
structive way. The direct coupling goes both in tune with the loop or 
against it. When 72 = ^'l = 2g the transfer is optimal. 



incoherent way, basically disrupting the Rabi oscillations, as 
we observed at the end of Sec. IVIBI G < g. 

Exceptionally, if the coherence transfer induced by pump 
and decay happens in a constructive, rather than independent, 
way as compared to the one induced by the direct coupling, 
Rabi oscillations and dressed states can be enhanced. This 
is the case of the optimally pumped configurations studied in 
Sec. I VI A| Let us analyze, for instance, the best situation, 
72 = Pi = and 71=^2 = 7> where G ~ V^g- In this case, the 
coherence transfer in the system can take place in a loop that 
involves not only the direct coupling but also incoherent pro- 
cesses, as sketched in Fig. [13] This results in the appearance 
of four different resonances in the system, that can be written 
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in terms of two pairs of new intermediate dressed states, 

|/±)^Z/^±5i|l,0)+/|0,l), (63a) 

|0±)ocZ/i±5i|l,0)+/|0,l), (63b) 
g 

named after inner and outer eigenenergies. They correspond 
to 7±/zi and 7±;z2, respectively, where zi.2 are those in 
Eq. (l47T i. The level scheme that these new coherence loops 
produce is plotted in Fig.|2d). In SSC, the dressed states ap- 
pear in the eigentransitions as: 

(Ta) = (|/_)(1,1|+/|0,0)(/+|), (64a) 

(Tb) = (|O_)(1,1|-/|0,0)(O+|), (64b) 

(Tc) = (|/+)(1,1|+/|0,0)(/_|), (64c) 

(To> = (|0+)(1,1 |-/|0,0)(O_ I). (64d) 

The normalization of these dressed states also depends on the 
coupling region. In SSC, inner stat es are norm alized by divid- 
ing by g{2g — 7) and outer by g{2g + 7), instead of ^/ig 
as was the case of the 4LS. These states do not have a bal- 
anced contribution from the bare states due to the saturation 
brought by pump and decay. The inner states have a smaller 
contribution from the first qubit than the outer states: 

l(1.0|/±>P = ^(l-^), (65a) 

|(l,0|O±)|2 = i(l + -l-). (65b) 

The 4LS dressed states correspond to an intermediate situation 
between the inner and outer dressed states in terms of popu- 
lations (equal to 1/2 for both qubit) and energy (see Fig.fTOl). 
This appears clearly in Fig. [T2]where the first qubit contribu- 
tion to each dressed state is plotted as a function of y/g. 

The two transitions composing each of the eigentransitions 
in Eq. ( |63] | have the same energy, which does not coincide 
with the average Hamiltonian value as soon as decoherence 
appears: {I±\H\ I±) = ±5R(zi)/[l - y/{2g)] and {0±\H\ 

0±)=±5R(z2)/[l + 7/(2§)]- 

When entering WC, the splittings between each pair of 
dressed states close, while the eigentransitions are linked each 
to a single state: 

(T^) = (|/_)(l,l |+/|0,0)(/_ I), (66a) 

(Tb) = (|0_)(1,1 |+/|0,0)(O_ I), (66b) 

(Tc) = (|/+>(1,1|+/|0,0)(/+|), (66c) 

(To) = (|O+)(l,l|+/|0,0)(O+|). (66d) 

Note that the inner peaks are always the most intense ones 
and that, when entering the WC, it is the "— " states that domi- 
nate the spectrum due to the saturation of the system into their 
limiting state | 1,0). 

Vin. CONCLUSIONS 

I described fully analytically the weak and strong coupling 
of two qubits in presence of dissipation and pumping. This 



configuration is realized by two two-level systems that com- 
mute. This was contrasted with other models of quantum cou- 
pling between two modes, namely, two bosons (that commute) 
and two fermions (that anticommute). 1 also addressed the 
case of a four-level system (4LS), which has the same level 
structure. 

Whereas the coupling of bosons, fermions and a 4LS mani- 
fest a conventional phenomenology, well known from the nor- 
mal coupling of two modes and differing only in bosonic or 
fermionic effective parameters, 1 showed that cross-Lindblad 
terms in the two coupled qubits give rise to a rich and com- 
plex new nomenclature of coupling. Namely, beyond stan- 
dard strong coupling (SC) where two dressed states emerge 
from quantum superposition of the bare states, I found that 
four dressed states, | I±) and | 0±), can be realized in the 
system with two different Rabi splittings, as a result of pump- 
ing and decay establishing new paths of coherence flow in the 
system. This is a new paradigm of strong coupling that leads 
to new regimes that 1 called "Second order Strong Coupling" 
(SSC), when the four dressed states are split (in two pairs), and 
"Mixed Coupling" (MC), when one pair of dressed states, the 
inner | /±), have closed. In this wider picture, 1 called "First 
order Strong Coupling" (FSC) the conventional type of strong 
coupling with only one pair of dressed states. These results 
show and spell out the considerable complexity of coherence 
transfer in other systems such as the Jaynes-Cumming model 
in presence of incoherent pumping, that are not amenable to 
analytical treatment. 
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Appendix A: Interactions and dephasing 

If the two qubits are close enough in space, they may 
interact. The energy levels are affected by the interaction 
(Coulomb interaction in the case of atoms or excitons) and 
the Hamiltonian must be corrected with the term: 

H^ = -X(yl<y\<yl(y2, (ai) 

where % is the binding energy between two excitations {x>Q 
corresponds to attraction). If J ^ O1.2, interactions have a 
negligible effect on the populations and statistics, but not on 
the spectral features, starting with Eq. ( l4b] i becoming ©1 1 = 
©1 + fi>2 — The dynamics in the presence of interactions 
is derived in appendix iB] although the expressions for spec- 
trum and dressed states are not presented in the text as they 
are not analytical in the most general case, and bring little to 
the discussion. The main effect that it produces is that upper 
and lower transitions can be identified thanks to the different 
energy of the upper transitions and this was used to check the 
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consistency of the results. In FSC and WC, the difference is 
already present due to the broadenings associated to each tran- 
sition and interactions only accentuate it. However, in SSC 
and MC, interactions break the symmetry between upper and 
lower transitions (e.g., cOa was linked to both upper or lower 
transition from inner states in Fig.|2d)) depending on the sign 

• if > 0, COa.d correspond to the upper transitions and 
(Ob.c to the lower transitions 

• if < 0, (Oa,d correspond to the lower transitions and 
COb.c to the upper transitions 

Naturally, strong interactions produce splittings between the 
peaks corresponding to upper and lower transitions, which 
may turn the spectral doublets into fully formed quadruplets. 

A second element that may be important in the system is 
pure dephasing. It may be caused, depending on the physi- 
cal system, by interaction with phonons or other external el- 
ements. This effect has proved to be relevant in state of the 
art semiconductor experiments Issll . Noting the rates of de- 
phasing of each qubit as 8i, 82, one must solve the following 
extended master equation from Eq. (|5]i: 



corresponding regression matrix M is: 



'^^ =^p+ ^ |(2c7/c7,pc7/cT,-c7/a,p-pa;cjO. (A2) 

!=1,2 ^ 



dt 



In the formalism, as a consequence, Eq. (iTal l becomes: 



(A3a) 



The dynamics in the presence of pure dephasing is again de- 
rived in appendix |B] and not discussed in the main text for 
the same reason as before, that it spoils analyticity of the re- 
sults and brings little to the discussion. The main effect of 
pure dephasing is to weaken the coherent coupling and is of a 
quantitative nature, not relevant enough to discuss it further in 
the present text. 



Appendix B: Quantum regression formula and time correlators 

In this appendix I obtain, thanks to the quantum regression 
formula, the first order differential equations for the one and 
two-time correlators required for the spectrum in Eq. ( fTOl ). 
The quantum regression theorem states that, once we find the 
set of operators C{^} and the regression matrix M that sat- 
isfy Tr(C{^}ifC) = L{A}^{T]A}Tr(C{A}0) for a general op- 
erator O, then the equations of motion for the two-time corre- 
lators (t > 0) read: 



^(o(f)q,}(f + T)) = Y.M{,,^{o{t)C{x}{t + T)) 

{X} 



(Bl) 



The most general set of operators for the two coupled qubits 



is C, 



{m,n,/j,v} 



Mmnfiv = i(0\ [m — n) + iOi2{ll 

iimllV 


-V)- 


iXimn 


-nv) 


(B2a) 
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- v)2 


M mnfiv = 
\—m,\—n,fiv 


P\mn — ix{l - 


ot)(1 - 
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-V), 


(B2b) 


M rnnjiv = 
mn.l — H.l — V 


Pi^V-ixil - 
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- v){m 


-«), 


(B2c) 


M mnfiv = 
m,\—n,\ — fi.v 


--2ig{v-ii){\ 
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(B2g) 



and zero everywhere else. 

We concentrate on computing g[^ \ as Gj'' can be obtained 
from it by exchanging the indexes 1 2. It corresponds to set- 
ting O = and having {m,n,iJ.,v} = {0, 1,0,0} in Eq. dBTT i. 
Other two-time correlators are linked to the one of interest by 
the regression matrix, and therefore are also needed to com- 
pute the spectrum. All correlators can be grouped in mani- 
folds that, following Refs H [H [M], I denote where 
k is the minimum number of particles that allows the corre- 
lators to take a nonzero value. As we can see schematically 
in the left part of Fig. [14] the first manifold ^ is composed 
of the correlators with {0, 1,0,0} and {0,0,0, 1}, linked only 
by the coherent coupling (red arrows). In the linear regime 
(very low pump), it is enough to solve the equations trun- 
cating in this manifold, which is the same as in the linear 

model and includes g[^\t). The next (and last) manifold 
where the two models differ, is composed of the correlators 
with {1,1,0,1} and {0, 1,1,1}. The links between manifolds 
are also of an incoherent nature, due to the pump (green and 
blue arrows). 

Gathering the four correlators in a vector 



v(f,f + t) 



/ (cT/(f)(T,(r + T)) \ 
(cT/(r)(T2(f + T)) 
((j/(f)cT/(7i(72(f + T)) 



(B3) 



their equations of motion read in matricial form: 
^v(r,f + T) = -Miv(f,f + T), 

with 



Ml = 



(B4) 



(B5a) 



CTj', with m, n, jj., V e {0, 1}. The 
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with 



FIG. 14: (Color online) Chain of correlators — indexed by {tj} = 
{m,n,il,v) — linked by the Hamiltonian dynamics with pump and 
decay for two coupled qubits. On the left (resp., right), the 
set IJi:=^ (resp., Uyt'-^) involved in the equations of the two-time 
(resp., single-time) correlators. In green are shown the first man- 
ifolds c/li and ^ that correspond to the linear model jT3], and 
in blue, the second manifold ^ and The equation of mo- 

tion (CTj' (0C{^} [t + t)) with rj e requires for its initial value the 
correlator (C{fj}) with {fj} G ,A^k defined from {tj} = (m,n,fl,v) 
by {fj} = {m+ l,n,il,v), as seen on the diagram. The thick red ar- 
rows indicate which elements are linked by the coherent (SC) dynam- 
ics, through the coupling strength g, while the green/blue thin arrows 
show the connections due to the incoherent quantum dot pumpings. 
The sense of the arrows indicates which element is "calling" which in 
its equations. The self-coupling of each node to itself is not shown. 
This is where COi 2 and Fi 2 enter. These links are obtained from 
the rules in Eqs. l lB2b . that result in the matrices Mi and Mq. The 
number of correlators needed to compute the spectrum is truncated 
naturally (with four elements) thanks to the saturation of the qubit. 



The general solution (for T > 0), v(f,f + t) = e'^^'^\{t ,t), 
leads to the correlator in Eq. (fT2l i with only four contributions, 
indexedby p =A, B, C andD. The frequencies (Op and rates Yp 
are the imaginary and real parts of the eigenvalues of the ma- 
trix M]. The dimensionless coefficients Lp{t) and Kp{t), also 
real, are linear combinations of the single-time average quan- 
tities, 
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(B8a) 



(B8b) 



This equation is similar to that of the linear model, only 
changing the effective fermionic rates into bosonic ones F, 
f/ [14]. On the other hand, the correlator with {1,1,1,1} 
(that is, (aj^aiaj <^2)(0)' which is the only one that is not 
zero in finds its expression separately, only in terms of 
the operators labeled {1,1,0,0} and {0,0,1,1}, through the 
equation 



d 

dt ' 



a/(Ti(Tja2)(f) 



-4r+((j/(Ji(j2^(j2)(f) 



(B9) 



In order to insert these averages in the expression ( fTOl ) for the 
spectrum, they must be either time integrated, in the sponta- 
neous emission case (SE), to give v^^ = [Q\{t,t)dt (and the 
coefficients L^^ and Kp^) or computed directly in the steady 
state (SS) to give v^^ = lim,^ooV(f ,f) (and the coefficients 
Lf, Kf). 



1. Two fermions 

The regression matrix Mi (cf Eq. iB5i ) truncates to its first 
submatrix. 



v(f,f) 



(cT/cT2)(f) 
((7((7/c7i(72>(0 
V(C7/(71(7](72>(0/ 



(B6) 



M 



1 — 




(BIO) 



recovering the exact mapping to the linear model, as discussed 
in the text. 



These mean values can also be found through the quantum re- 
gression formula, setting the operator 0=1 and deriving the 
corresponding regression matrix from the rules in Eq. ( IB2l i. 
They can also be grouped in manifolds, denoted All 
nonzero single-time average quantities, all of them needed to 
compute the spectrum, are shown in the right part of Fig. [14] 
Note that ((7j^(7|(7i(72)(0 — 0. The correlators in ^ can be 
obtained independently solving a second matricial equation: 



du{t) 
dt 



-Mou(f)+p 



(B7) 



2. Matrix of coherences 



The matrix of coherence in Eq. ( |541 ) can be found by rewrit- 
ing the regression matrix Mi in the basis of projectors for 
each transition, presented in Eq. ( |53] i: Mcoh ~ A^'MiA with 
A some transformation matrix. For this reason, they have 
the same eigenvalues, which are the resonances of the sys- 
tems. Since the eigenstates of Meoh are the transitions be- 
tween dressed states, their study enables us to reconstruct 
them even in a highly dissipative environment. 
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Appendix C: Cross Lindblad terms 

In the case where the four levels of the system do not cor- 
respond to two qubits but are the levels of a single entity ll32ll . 
the corresponding Liouvillian of the system is most conve- 
niently written in terms of the four-level operators in Eq. (|6]l. 
The final master equation in that case can be found from 
in Eq. (|5]l by removing cross terms that appear in the two cou- 
pled qubits; 



dp 
dt 



- ^p-Yi ihpu] + uipl\) - Yiikpul + uipll) (CI) 



-Pl{u\ph+l[pUi)~P2{u2pl2 + llpU2) (C2) 

-5\{l[l\pu\u\ + u\uipl[l\) 
ih ^.,t., , .,t. 



-d2(l2hpU2U2 + U2U2pl2h) ■ 



(C3) 
(C4) 



Such cross terms represent effective couplings between the 
coherences, in the case of pump and decay, and the states, 
in the case of pure dephasing, associated to each qubit. Ev- 
idently, they can only have consequences on the dynamics 
when the two qubits are also directly coupled. The rules to 
add to those in Eq. ( IB2b in order to obtain a new general re- 
gression matrix M + 8M are: 

dMnmtiv ^P\{m + n- - v)^ +P2(M + V- l)(m-m)^, 

mnllV 

(C5a) 

8M mniiv =-(jU-v)^[(Pi + 5i)m«-ri(l-m)(l-«)], 

(C5b) 

5M mniiv =-(m-n)^[(P2 + 52)AiV-f2(l-jU)(l-v)]. 

(C5c) 

With these new rules, the matrix for the one-time correlators, 
Mo, is still the same as in Eq. ( |B8bt and, therefore, the mean 
values in u(f), Eq. ( IB8al ). do not differ from those for two 
qubits. The two-time regression matrix Mi changes into Mi + 
5Mi with 



5Mi 



Appendix D: Dressed states in the manifold picture 

We can give an alternative derivation and intuitive interpre- 
tation of the results in section [VAl bv assuming complex fre- 
quencies Oi 2 — > iii,2 = fi^i,2 ^ '71.2/2 in Eq. (O. This leads 
to the complex dressed frequencies: 
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From here, we can compute positions and broadenings of the 
four possible contributions to the spectra. The positions are 



given by the difference in energy between the levels involved 
in the transitions while the broadenings are given by the sum 
of the imaginary parts (the uncertainties). If we apply this 
principle between the levels in Eq. dPlI ). we find the same 
result than in Eq. (I34li-(l35]l. This allows us to identify the 
contributions A, B, C, D to the SC spectrum with the mani- 
fold transitions in the order that we now explain. In Fig.|2b) 
we can follow the four possible transitions in the manifold 
picture. The lower transitions (in blue), from | ±), coincide, 
respectively, with the expressions for peaks A and D, given by 
Eq. ( |34] i. The upper transitions (in red), towards | ±), coincide 
with the expressions for peaks B and C, given by Eq. ( l34l i. The 
upper transitions have a larger broadening than the lower due 
to the addition of the uncertainties in energy of the levels in- 
volved, brought by the spontaneous decay. 

In order to understand better some features of the spectra of 
the steady state under incoherent pump in the different regions 
that we have defined in section IV Bl we will now attempt to 
push the manifold method — adequate for vanishing pump — a 
bit further, to the case of non-negligible pump. If we combine 
their effects in the non-hermitian Hamiltonian as 



P1+P2 



n± = (Oi-ir+±R 
71 + 72 



ITD 



2o)i 



(D2a) 
(D2b) 
(D2c) 



and apply their sum and subtraction to obtain positions and 
broadenings of the spectral contributions of each of the four 
transitions: 
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(D3a) 
(D3b) 
(D3c) 
(D3d) 



we do not obtain the right results of four — in principle 
different — peaks. This discussion brings us back to the naive 
association of SC with 7^ we discarded in Sec.HIIl 

It is evident that using only the single-time dynamics fre- 
quency R^^^ to describe the splitting between the dressed 
states, cannot give the variety of situations than the two com- 
plex parameters zi 2 give. 

In the case of a 4LS, the results obtained in Eq. dJTT i are 
comparable to those in Eq. ( |D3t , even though the positions 
of the peaks are given by R instead of /?'™. The broaden- 
ings are exactly the same. We can conclude that without the 
cross Lindblad terms, the coupling acquires a bosonic charac- 
ter: the correlations between levels are relaxed from those of 
two qubits into those of one (larger) four-level system where 
the manifold picture holds. 
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Appendix E: Second order correlations functions 



In terms of the following vector of mean fluctuations. 



In this appendix, I obtain the second order direct and cross 
correlation functions defined, respectively, as 

Cf '(f ,r + T) = (cj/(r)cj/(r + t)cji (r + T)a, (r)) , (Ela) 
Gf^{t,t + t) = (c7/(f )c72^(r + T)c72(r + t)(7i (f)> . (Elb) 
These two quantities are proportional to the probabilities of 

(2) 

two T-delayed emissions from the same qubit {G\ ) or two 

(2) 

different qubit emissions (GJj )• 
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FIG. 15: (Color online) Steady state values of (solid blue) 

(t) (dashed purple) at resonance, (a) The system is in FSC: 71 = 
g, 'Y2 = g/2. Pi = g/2 and P2 = g/\Q. (b) In MC with optimum 
coupling G = \/2g: Ji = f 2 = and P\= J2 = 2g (see Fig. [Hi. 
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(E5) 



the steady state normalized correlations functions (for T > 0) 
read 
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Gf2^(0,T)_^ [g-Mo^fSS]2 



' 12 
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e[0,l], (E6a) 
e [0, 1] (E6b) 



where [x]i means that we take the first element of the vec- 
tor X. At infinite delays, the two functions reach the value of 
uncorrelated emissions. 



(E7) 



In order to compute such two-time correlators, we apply the 
quantum regression formula presented in an alternative form 
to Eq. dHB: 



{X} 



(E2) 

where 2 is a second arbitrary operator The equations we must 
solve setting O = and Q~ Gi, are the same as those used to 
compute the mean values in Eq. (IB7b but now with the vector 
of correlators 



w(f,f- 



-T = 



/(C7|(0CT/C7l(r 
(j/(0cjlcJ2(f- 

a/(r)c7/a2(f ^ 
V(cr/(f)c7iC72^(f- 



-^)cTi(0)\ 

-T)(Ti(f)) 
-T)(Ti(f)> 

-T)cTi(O)/ 



and the vector ni (f )p instead of p. In the steady state (f = 0), 
the solution for w(0, t) requires knowing the initial condition: 



w^^ = w(0,0) = 




(E4) 



At zero delay, it is evident that g'''(0) = 0, as two excita- 
tions cannot exist in the same qubit. At intermediate delays, 
the emission presents antibunching: g'^) (q) <g''^\x), the sec- 
ond excitation from the same qubit has more probability to be 
emitted after some delay (see Fig. [TSl l. At zero delay. 



8>) 



n\n2 



(E8) 



It always is possible a second emission coming from a differ- 
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(E3) ent qubit (the minimum value for g\2{,l = 0) is 1 /5). When 

(2) 

the qubit behave as independently, g j2 (t) = 1 all along. 

Being the second order correlation functions given by basi- 
cally the same equation as the populations and depending only 

f21 f2) 

on them, g\ (t) and ^12 ('f') are the same for the two coupled 
fermions, qubits and 4LS. It is not possible to distinguish from 
their dynamics the four coupling regions described in Table U 
appearing in the case of coupled qubits. 
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